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Numerical results for the polaron dispersion are presented for an arbitrary number of space di- 
mensions. Upper and lower bounds are calculated for the dispersion curves. They are rather close 
to each other in the cases of small electron-phonon couplings usual for real polar materials. To 
describe the dispersion in other materials, we suggest a simple fitting formula which can be applied 
at intermediate values of the Frohlich electron-phonon coupling constant. Its validity is approved 
by the comparison with direct calculations and previously obtained results. This makes our results 
not only reliable and highly accurate but also easy reproducible. 
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I. INTRODUCTION 

The electron-phonon interaction influences the proper- 
ties of charge carriers (electrons) in polar semiconductors 
or ionic crystals. An electron polarizes a medium and, 
being surrounded by a cloud of virtual phonons, is cap- 
tured by a self-induced potential which can move in a 
material. Such a quasiparticlc is called a polaron. The 
larger the value of the Frohlich electron-phonon coupling 
constant a, the more pronounced are polaron effects. In 
particular, the electron-phonon interaction results in an 
electron binding energy, in a renormalization of its mass 
and in a nonparabolic energy-momentum dependence. 

In the present paper, we study the large polaron disper- 
sion law that is the dependence of the polaron ground- 
state energy E{a, Q) on the total polaron momentum 
Q. The technological progress in man-made structures 
has caused a rapidly increasing literature on systems of 
reduced dimensionality, in particular, on the effects of 
electron-phonon interaction in quantum wells, wires and 
dots. Thus, the dimensionality D of space in our study 
may be different: I? = 1, 2, 3. 

Here and in what follows all quantities are dimension- 
le ss, the e nergy, mass and length units being Uuj, m and 
y^h/2muj, respectively, where ui is the longitudal optical 
(LO) phonon frequency and m is the free electron band 
mass. For instance, the polaron kinetic energy /2mpoi 
at small momentum P is written down as /m(a) in 
our dimensionless units where is the polaron effec- 
tive mass and m{a) = rripoi/m. 

The literature on polaron is enormous but it concerns 
mostly the energy E{a,0) of the bulk {D ~ 3) polaron 
at rest. Nevertheless, the first results were obtained on a 
bulk polaron dispersion law even in a very early paper by 
Frohlich, Pelzer and Zienai>i where the authors used the 
first order of the Brillouin-Wigner perturbation theory. 
Later Whitfield and Puff^ suggested an improved version 
of the polaron energy-momentum relation in the weak- 
coupling regime. For these and other early papers see also 



the review article by Appel^ where qualitative considera- 
tions on the behavior of the dispersion curve were given. 
The results obtained demonstrate that the bulk polaron 
energy-momentum relation is quadratic for small Q but 
then bends over and becomes horizontal when the energy 
approaches the continuum edge Ec which is reached at 
some finite value Qc of the polaron momentum. At this 
momentum the moving polaron energy Ec = E{a, Qc) 
exceeds the ground-state energy E{a, 0) of the polaron 
at rest exactly by the energy of a free phonon Hlu (which 
is just unity in our notation). 

Below the continuum edge Ec = E{a, 0) + 1 (that is at 
Q < Qc) the ground-state energy E{a, Q) is an isolated 
and well defined eigenvalue. There are some important 
rigorous results concerning the properties of the disper- 
sion E{a, Q) (see Ref.^>^) 

• E{a, Q) is a real analytic function of a and Q for 
0<a<oo,Q^<l. The former constraint on Q 
can be removed totally for D = 1,2. For D = 3 the 
domain of Q can be extended up to a finite value 
Qc, where the energy reaches the continuum edge. 

• E{a, Q) decreases with a and increases with Q be- 
low the continuum edge. 

• The inequality E{a,0) < E{a,Q) holds for g 7^ 
and < a < 00. 

• The upper bound is given by the inequality 

E{a, Q) < mm[E{a, 0) + Q^ E{a, 0) + l]. (1.1) 

The simplest and seemingly the most natural way to 
describe the dispersion is the first order of the Raleigh- 
Schrodinger perturbation theory (RSPT). For D = 1,2,3 
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it results in the well-known formulae: 

TT 



E{a,Q) = Q'~a- 



for D^l, (1.2a) 



E{a,Q) ^ ~ aK{Q^) for D ^2, (1.2b) 
E{a, Q) ^ ^arcsmQ ^ ^ ^ ^c) 

Here K{m) is the complete elliptic integral of the first 
kind. 

The results of the RSPT contradict the rigorous prop- 
erties of the dispersion. For instance, all three functions 
of Eq. p.2p have maxima at some momenta Qm < 1 so 
they are decreasing functions in the region Qm < Q < 1. 
Moreover, the expressions for = 1,2 diverge at Q = 1, 
so the RSPT evidently fails to work near this value. At 
small Q <C 1 RSPT leads however to correct (to the first 
order in a) parabolic functions of the effective mass ap- 
proximation. 

After the above cited paper— different variational ap- 
proaches were developed to calculate the polaron disper- 
sion. The paper by Lee, Low and Pines^ should be men- 
tioned among earlier articles on polarons. For sufficiently 
large Q their variational result is weaker than the nonan- 
alytical upper bound (jl.ip , the dispersion curve intersects 
the continuum edge and the deviation from the correct 
result is of the order , which does not provide the nec- 
essary accuracy of calculations. 

More advanced variational calculations were performed 
by Larseni, and by Warmenbol, Peeters and Devreese^i^. 
These and some other papers will be discussed in more 
detail later. 

Naturally, only upper bounds could be obtained with 
the variational methods of these papers. It is much more 
difficult to derive a lower bound, and it was found by 
Lieb and YamazakiiS but only for Q = 0. The general 
procedure to obtain the Q— dependent lower bound was 
developed by Gerlach and Kalinali. 

In the present paper, we combine the variational upper 
bound obtained with the expansion of the trial wave func- 
tion in numbers of virtual phonons with the lower bound 
of Ref.^^. In principle, our variational upper bound could 
lead to exact solutions but in practice we have to cut the 
expansion and work with an approximation obtained this 
way. For small values of the electron-phonon coupling 
constant a which are common for most of the polar ma- 
terials the corridor between these two estimates is very 
narrow, so that we can pretend to finding numerically 
exact solutions. 

The lower bound gives too poor results for intermedi- 
ate values a ~ 1. Besides the huge numerical job does 
not allow us to reach the necessary accuracy with the 
upper bound. To overcome these difficulties, we suggest 
simple fitting formulas to calculate the polaron disper- 
sions in different dimensions for intermediate values of 
the coupling constant. This makes our results easy repro- 
ducible and reliable, which is demonstrated while com- 
paring them with these by other authors. 



II. BASIC EQUATIONS 

The starting point is a Hamiltonian of Frohlich type 
in the form which makes use of translation invariance 
to perform a projection onto a subspace of fixed total 
polaron momentum: 



HiQ) : 




(2.1) 



The electron coordinates were eliminated with the well- 
known Lee-Low-Pines canonical transformation^. This 
reflects the conservation of the total polaron momen- 
tum Q which is a c-number in Eq. (j2.ip . Here k, gk,ak, 
and aj^ are the wave vector, coupling function, and the 
annihilation- and creation operator of the phonon under 
consideration. The coupling function is defined as 
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Vk^- 



(2.2) 



where V is the quantization volume and is a number. 
For D = 3 and ^3 = 4 or _D = 2 and ^2 = 2 one recovers 
well-known models (see, e.g. RefJ^). At first glance, the 
case D = 1 has to be excepted as, according to Ref.— , 
^1 diverges. Nevertheless, the coupling is physically 
interesting for D — 1 and finite ^1 - either in the sense 
of a regularized version of a polaron model, as discussed 
in Ref li^, or in the sense of an effective model within the 
theory of the bulk {D = 3) free polaron model-^-. We 
choose ^1 = 1 without loss of generality. 

Summations in final formulae are replaced by integra- 
tions following the conventional rule: 
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.9k I' Fik) 



-F(k). 



(2.3) 



The scheme to calculate the upper bound for E{a, Q) 
was presented in Refi^. The variational principle of 
Ritz was used. Choosing an adjustable, normalized 
wave function | $ > and calculating the minimum of 
< $ I i?(Q) I $ >, one gets an upper bound z = z(a, Q) 
to E{a, Q). Our trial function is of the type 



$ >:= C 



0> - JZ^kSkfll, I 0> 



(2.4) 



y^gkffk'-Bk.k'Q^Q^/ I > 



k,k' 



- E 5kffk'5'k"^k,k',k"akak/a[// | > 

k,k',k" 



the two- and three-phonon amplitudes being totally sym- 
metrical functions. The mean value < <f> | -ff (Q) | $ > 
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can be readily calculated from which we obtain the min- 
imizing equations for the amplitudes. The more phonon 
amplitudes are included, the better is the upper bound. 
In principle, one can include an arbitrary number of 
phonons arriving at a subsequently large number of equa- 
tions for the corresponding amplitudes. 

In practice, one has to cut expansion (|2.4p . For exam- 
ple, restricting ourselves to the one-phonon amplitude we 
arrive at the Brillouin-Wigner perturbation result in the 
first order in a: 



z = u- - 



E 



I5k| 



(Q-k)2 + l-z' 



(2.5) 



If we take into account the two-phonon contribution, 
the system of subsequent equations takes the form: 



(2.6) 



and 



[(Q - k)2 + 1 - z] = 1 + 2^ |.9k'|'5k.k', (2.7) 



2 [(Q - k - k')^ + 2-z] BkM' = Sk + Sk' . (2.. 



The quantity i?k.k' can be readily found from Eq 
Inserting it into Eq. (|2.7[) and defining the function 



<x;z) = 



l5k 



f (x - k')2 + 2~z 

we arrive at the equation for the quantity B]^: 

[(Q - k)2 + 1 - z - d{Q - k; z)] Bk 
IffkfSk' 



(2.9) 



1 + E 



f (Q-k-k')2-|-2-z" 



(2.10) 



To keep the necessary accuracy of the solution, we it- 
erate it twice and truncate the series: 
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iV(Q-k;z) 



(2.11) 
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f N{Q - k; z)[(Q - k - k')2 + 2 - z] ■ 



where 



7V(x; z) ^ + 1- z - d(x; z). 



(2.12) 



Inserting Eq. (|2.12p into Eq. (|2.6p we obtain the upper 
bound: 
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iV(Q-k;z) 



(2.13) 
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Lgkhs: 



2|„ |2 



k,k 



7V(Q - k; z)[(Q - k - k')2 + 2- z]N{Q - k'; z) ' 



Recall that the dispersion curve starts at the minimal 
value E{a,Q) and reaches its maximal value at the con- 
tinuum edge. The exact expression for the latter is as 
follows: 

E^ = E{a,Q) + l. (2.14) 
In our variational estimates this formula is modified: 



E, 



c.n ph 



E, 



(n-l) ph 



(a,0) + l, 



(2.15) 



where the lower indices show the number of phonon 
amplitudes taken into account. Here n > 1 and 
£^Oph(Q^,0) = (no electron-phonon interaction). The 
interpretation of Eq. (|2.15p is clear: as we work with a 
fixed number of phonon amplitudes and one phonon be- 
comes free at the continuum edge, the number of phonons 
which contribute to the polaron state decreases by one. 

Therefore, in the one-phonon approximation we obtain 
for the continuum edge a rather trivial and poor estimate 
Ec,ip\i = 1. If a three-phonon amplitude is included, 
we describe z{a,Q) with high accuracy (e.g., correct up 
to terms of order near the continuum edge and up 
to terms of order for the small-Q behaviour). This 
would be enough for reasonable estimates in the weak 
coupling regime. The only problem is that the numerical 
job becomes time consuming and could be done only for 
the case D = 1 (see Refj^). 

To conclude this section we mention a qualitatively dif- 
ferent behavior of the dispersion curves for Z) = 1, 2 and 
D — i. The continuum edge is an asymptote for D = 1,2, 
and is approximated from below as Q — > oo, whereas 
for D = 3 the dispersion does meet the edge at a finite 
value of Q = Qc, which depends on a (see Refsi^ i^^'^^ ). 
This distinction between the low-dimensional and bulk 
polaron dispersions is explained by properties of the in- 
teraction potential and may be readily understood. The 
total polaron energy z is presented in the r.h.s. of Eq. 
(|2.5p as a sum of a positive free polaron kinetic energy 
and a negative interaction energy proportional to a. 
At D = 1,2 the latter becomes infinitely large when the 
total energy z approaches the continuum edge (unity in 
this approximation) because the integral in Eq. (|2.5p di- 
verges. To compensate this and to keep the total energy 
finite, the kinetic energy (and the polaron momentum) 
should tend to infinity as well. This is not the case in 3D 
where the interaction energy stays finite at the contin- 
uum edge and so does the limiting value Qc of the total 
polaron momentum. 



III. THE ID CASE 

With only one-phonon exchange taken into account 
formulae (|2.2p and (|2.5p result in the following equation 
for the upper-bound z{a, Q) at Z? = 1: 
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2^/l~~ 



(3.1) 
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If we put z = in the r.h.s. of Eq. (|3.ip . we immediately 
arrive at the expression for the first order of RSPT given 
by Eq. ^L2ah . 

An upper bound obtained with the inclusion of the 
two-phonon amplitude can also be presented in a closed 
analytical form: 
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FIG. 1: Dispersion of ID polaron in GaAs 

General formulae and some numerical results for a 
three-phonon exchange as well as a lower bound can be 
found in Ref.^^. Here they are presented in Fig. [T] for 
a very small value of the coupling constant a = 0.068 
corresponding to a real material (GaAs). 

IV. FITTING FORMULA 

To give an idea of typical values of the Frohlich cou- 
pling constant for other materials, we present here some 
experimental data from Ref.^: a = 0.022 for InSb, 
a = 0.123 for AlAs, a = 0.65 for ZnS, a = 1.84 
for AgCl. It is well known that intermediate values 
a 1 and even larger can be tackled within conven- 
tional Raleigh-Schrodinger perturbations. This becomes 
possible because the coefficients e„ in the expansion 
E{a,0) = ^n^nCt" of the polaron ground-state energy 



decrease very fast with n so that the expansion is per- 
formed, roughly speaking, in powers of a/10 or so. One 
can notice this tendency (although it is not proved rig- 
orously) looking at Eq. (|6.7p . The famous Feynman 
approximation^* clearly demonstrates this property as it 
follows from our calculationsi^ of first twelve (sic!) coeffi- 
cients of the weak-coupling expansion for the 3D-polaron. 
Besides, the radius of convergence Ra of the RSPT for 
the bulk polaron ground-state energy was estimated in 
a number of papers: Larsen^S found Ra ~ 6.2 6.5, 
Klochikhin^i argued that Ra ^ 3.4, our estimatei^ gave 
Ra 6.9. 

Anyway, this explains why the weak coupling expan- 
sion can be applied for intermediate values a ^ 1 or even 
larger. On the other hand, our lower bound gives very 
poor results for the intermediate coupling. In addition, 
the two-phonon approximation works not so well at these 
values of the coupling constant a. The reason is that the 
edge point i?(Q;, 0) -I- 1 is then calculated within the one- 
phonon approximation, which is certainly not enough to 
reach appropriate accuracy. 

Taking as an example ID-polaron, Eiph (a, 0) is defined 
as z of Eq. at Q = 0: 

Eiph{a,0)^l - Eiph{a,0) = -a^. (4.1) 

Expanding Eiph{a,0) in powers of 7ra/2 we obtain 
Eiph{a,0) = -{aTT/2) + {l/2){aTT/2)^+0{af. Here the 
second order coefficient is positive and equals one-half, 
while in the correct resuli^i^ 

i^Wa,0)=-(f)-(^-l) {^)\0iar 

it is small (-0.06) and negative. From the point of view of 
diagrammatic technique for polarons^^ it means that only 
the disconnected Feynman diagram is taken into account 
in the one-phonon approximation. In this approximation 
the contribution of two connected diagrams is skipped 
which otherwise would compensate rather a large posi- 
tive coefficient to give a small negative residue -0.06. An 
analogous fine tuning happens for flat {D = 2) and bulk 
{D = 3) polarons as well. 

This explains why one needs to include three-phonon 
amplitudes to describe not only weak couplings but also 
intermediate values a ^ 1. However, as was mentioned 
earlier, three-phonon calculations can be performed only 
in D = 1, and the numerical job becomes enormous in 
other spatial dimensions. So we propose now a simple 
formula to fit the exact dispersion curves. We demon- 
strate its validity for Z? 1 by the comparison with the 
three-phonon approximation. In other sections of the 
current paper the fitting formula will be used instead of 
the absent numerical three-phonon calculations. 

Let us stress that we are not going to proceed to the 
strong-coupling limit. We are still dealing with the weak 
couplings and our goal is to restore what is lost in ap- 
proximations with partial summation of the perturbation 
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series, that is to regain a possibility to use the weak- 
couphng resuhs for the intermediate values of a. 

We start with the first term of the Brillouin-Wigner 
perturbation series (j3.ip which reproduces the general 
behavior of the dispersion curve. As we already know, 
its main disadvantage is the lack of accuracy, especially 
in calculating the value of the continuum edge. Thus, we 
propose to remedy this weak point introducing correction 
terms " by hand" . 

First of all, we replace the first order expression for 
the ground-state energy (which equals —an/ 2 in this 
case) by its exact value E{a,Q). Then we also replace 
the approximate edge point which is unity by its ex- 
act value 1 -I- E{a, 0) in the propagator of the r.h.s. of 
Eq. (I3.1|l . Note that actually this step is not an ap- 
proximation: the energy can be arbitrarily shifted in the 
denominators of the Brillouin-Wigner expansion, as was 
mentioned in Ref.^'^. This allows us to keep the correct 
gap between the zone's bottom and the edge point. Fi- 
nally, we scale the total momentum Q by the factor bi 
to obtain the correct effective mass behavior of the type 
E{a,Q) « E{a,0) + Q'^/m{a) at small Q. This way we 
arrive at our fitting formula: 



E{a,0) 



^l + E{a,0)-f' 



(4.2) 



m(a) 



1 - 



E{a,0) 



Later we apply a similar procedure to the cases D — 2,3. 

At small Q the correct parabolic behaviour is guaran- 
teed by the very construction of Eq. (|4.2|) . When Q is 
large it leads to the asymptotic behaviour 



f^l + E{a,0)- 



(4.3) 



Thus, the dispersion curve approaches the correct con- 
tinuum edge rather fast. 

The exact expressions for the polaron ground-state en- 
ergy E(a, 0) and its effective mass m(a) are unknown but 
they can be calculated numerically with the help of the 
perturbation series. The first two terms of the latter were 
found in Ref.^— and a few next terms were calculated by 
Khomyako\*2i: 



E{a,0) 



0.060660 f^V- 0.00844437 (^)' 



V 2 



-0.00151488 — 
2 



an 



m{a) = l + 0.5(^—j + 0.191942 (^—j 



-0.0691096 



V 2 



(4.4) 



The results for a = 0.5 are presented in Fig. [51 where 
numerically exact ground-state energy E(a, 0) and the 
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FIG. 2: Dispersion of ID polaron for a = 0.5 



edge-point 1 + E(a^ 0) are shown by solid thin straight 
lines. 

Our fitting formula is close to the so-called improved 
Brillouin-Wigner perturbation theory of Rcf.^"^ which ap- 
pears to be equivalent to the one-phonon Tamm-Dancoff 
approximation^. Previously, this approximation was 
used for polarons in = 2,3 (the discussion and refer- 
ences will be given later) . With these ideas being applied 
to the polaron in ID, we obtain the following expression: 



TD 



anjl 



\J\- an 12 - fTD 



(4.5) 



In comparison with our Eq. (|4.2|) only the first order in a 
is taken into account for the polaron ground state energy 
at rest and the effective mass (which gives bi = 1). The 
curve fro is also plotted in Fig.O 

Comparing with the three-phonon calculation one can 
notice that our fitting curve provides us with an excel- 
lent result which is also better than the Tamm-Dancoff 
approximation. This is not surprising because we used 
in our construction the RSPT expansions to rather high 
orders which work quite well for the intermediate values 
of a. Thus, analogous fitting formulas will help us in the 
next sections where it is not possible technically to per- 
form calculations within the three-phonon approximation 
to proceed to the region of intermediate couplings. 



V. THE 2D CASE 

With the one-phonon exchange taken into account (the 
same Eqs. (|2.2p and (|2.5p taken for D — 2) we arrive at 
the equation for the upper bound: 



K 



1 - z 



(5.1) 



If we put z = in the r.h.s. of Eq. (|5.ip . we get 
the expression for the first order of the RSPT given by 
Eq. (fOb]) . 



ergy and effective mass can be taken from Ref] ^^'^^ : 
E{a,0) = -|a- 0.063974a^ 



m 



(a) = 1 + -a + 0.127235a2. 



(5.4) 
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FIG. 3: Dispersion of 2D polaron in GaAs 
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The resuhs obtained with this formula and with the 
two-phonon exchange as well as the lower bound are pre- 
sented in Fig. [3] for a = 0.068. Numerical calculations 
for the three-phonon exchange are very complicated and, 
therefore, will not be performed in this case. However, 
for such a small value of a the gap between the upper 
and lower bounds is very narrow. 

To tackle intermediate values of a, we apply the same 
idea to construct a fitting formula which reads in 2D as 
follows: 



/ = 



xK 



E{a,0) 



{h2Q? 



hi 



{b2QY + l + E{a,Q)- f 
1 E{a,Q) 



(5.2) 



m{a) 



1 - 



4-£;(a,0) 



Here again the correct parabolic behaviour at small Q 
is reproduced by the very construction of the fitting for- 
mula. The asymptotic behavior of the complete elliptic 
integral K{m) ^ —(1/2) ln(l — m/16) when m ~* 1 leads 
to the following asymptotics of the dispersion curve at 
large Q: 



/wl + i;(a,0)-16(62Q)'exp[- 



\E{a,Q)Y 



(5.3) 



Thus, the dispersion curve approaches the continuum 
edge extremely fast. 

The expressions for the 2£'-polaron ground-state en- 



FIG. 4: Dispersion of 2D polaron for a = 0.5 

The Tamm-Dancoff approximation for the 2D polaron 
was calculated by the Antwerp group^. Again, the for- 
mula is very similar to our Eq. (j5.2p , where 62 is replaced 
by unity and E{a, 0) is put equal to — a7r/2: 



ho = Q'- 
xK 



+ aTT/2 - Jtd 



^Q^ + l-a7T/2-fTD 
The numerical results are shown in Fig. [J] 



(5.5) 



VI. THE 3D CASE 

With Eqs. (|2.2p and (|2.5p the result of the one-phonon 
approximation for the bulk polaron {D = 3) reads as 
follows: 



^2 

z = U — — arcsm ■ 



(6.1) 



This equation appeared at first in Ref.^. Putting z — 
in the r.h.s. of Eq. (|6.ip we get the expression for the 
first order of the RSPT given by Eq. pT^ . 

Being expanded in powers of the coupling constant this 
reproduces at Q = the exact results to the first order 



z(a,0) = -a + 0{a'^). 



(6.2) 



For the edge point we have again Zc — I. Inserting this 
value into Eq. (j6.ip we obtain the equation for the value 
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Qc at which the edge point is reached: 



(6.3) 



In this approximation the curve z{a, Q) approaches the 
edge point z{a, Qc) = 1 as an inverse parabola: 
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FIG. 5: Dispersion of 3D polaron in GaAs 

The curve for the upper bound which we calculated 
within the two-phonon approximation for GaAs is shown 
in Fig. [5l In the whole range of momentum Q it is al- 
most undistinguishable from the lower and upper bounds 
within the one-phonon approximation. Thus, it can be 
considered as a (numerically) exact result. All three 
curves are presented near the edge point in Fig. [S] The 
curvature of the dispersion near the edge point is very 
large at small a, as it follows from Eq. (|6.4p . and the in- 
verse parabola described by this equation can hardly be 
seen in these plots. 

Again the gap between the edge point and the zone's 
bottom is larger than the phonon energy (unity) in the 
one-phonon approximation. There is no physical reasons 
for that, and this disadvantage was removed in Refi^. 
However, the authors obtained wrong weak-coupling ex- 
pansions for the polaron energy and effective mass at 
small momenta. A modification of Eq. (j6.ip was given 
by Klochikhin^ who studied the 3D polaron dispersion 
in the scope of the perturbation theory and took into ac- 
count two-phonon amplitudes. This way he arrived at 
the following equation: 



TD 



— arcsm — , 



a 



(6.5) 



TD 



which coincides with the Tamm-Dancoff approximation^. 
This way he obtained the value Ec ~ \ — a for the con- 
tinuum edge in the first order in a which is reached at 



Qc = 1.16. The correct energy shift was the advantage 
of these calculations in comparison with previously made 
in Ref..^. 
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- Lower bound k{a,Q) 
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FIG. 6: Dispersion of 3D polaron in GaAs near the edge point 

For intermediate values of a (we present here exam- 
ples with a = 0.5, 1) the two-phonon exchange contri- 
bution does not provide us with impressive results. The 
three-phonon contribution requires a very large amount 
of computational work and, therefore, will not be given 
here. Besides, such values of a are outside the domain 
admissible for the lower bound, although the polaron is 
still in the weak-coupling regime. So we again construct 
our fitting formula which now takes the form: 



/ = {hQf + 
1 



hQ 



m(a) 



1 



— — — — arcsm — , , 
hQ ^{b^QY + l + E{a,{))~ f 

E{a,0) 



6-2E{a,0) 



(6.6) 

a this equation 



Note that at 63 = 1 and E{a,0) - 
reduces to Eq. (|6.5p . 

The bulk polaron ground state energy and its effective 
mass are given by the known perturbation series (the 
second order in a was found in Refi^ and the third order 
was calculated in Refj^^i^^: 

E{a,Q) = -a- 1.59196(a/10)2-0.80607(a/10)^ 
m{a) = l + a/6 + 2.36276(a/10)^ (6.7) 

The maximal value Qc of the polaron momentum can 
be found from the fitting formula (|6.6p if we put / = 
1 -|- E{a, 0). Then we arrive at the equation for Qc'. 



Qc = Qc/bs, 

Qc {q'c-l-E{a,0)) 



--Eia,0). (6.8) 



The fitting dispersion curve f{a,Q) approaches the 
edge point also as an inverse parabola at Q « Qc, Q < 
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f{a,Q) « l + E{a,0) 



{qc-q)\ (6.9) 



where q — b^Q. 

This way we found = 1.02 for a = 0.068, Qc = 1.11 
for a = 0.5, and Qc — 1.20 for a = 1. The fitting curves 
obtained are shown in Figs. [7] and [5] together with one- 
and two-phonon calculations. Being apphed at a = 0.068 
this procedure leads practically to the same results as 
direct calculations. 



0.4 
0.2 



Upper bound z(a,Q), 1 phonon 

Upper bound z(a,Q), 2 phonons 

- - Fitting curve 
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FIG. 7: Dispersion of 3D polaron for a — 0.5 

The results of the Monte-Carlo calculations^!^ are 
shown also in Fig. [S] The relation k — QV2 is used 
between our momentum Q and their momentum k (free 
electron energy reads as follows: E = Q^ = k'^/2). We 
see excellent agreement between our results almost in the 
whole range of momentum except for the vicinity of the 
edge point. The Monte-Carlo calculations give the value 
Qf^ w 1.285 while our fitting formula gives Q;?" = 1.20. 
This discrepancy of 8% is responsible for the deviation 
of these curves. The authors of Refsj ^^'^^ reported that 
their calculations were "numerically exact" and the er- 
ror bars were smaller than the size of points at the plot. 
Previouslyji^ we criticized these statements, in partic- 
ular, because the results obtained give the wrong coeffi- 
cient even in the second order in a for the polaron ground 
state energy weak-coupling expansion. As we could see 
above, the correct second order calculations are crucial 
for the adequate description of the dispersion near the 
edge point. 

The polaron dispersion /(a, Q) was also calculated by 
Larseiii who combined for his variational ansatz the one- 
phonon Tamm-Dancoff approximation and the Lee-Low- 
Pines transformation^: 



—— — — — = — F a,g), 
jL + a + Q^ n 



F(a,Q) 



1 /L + a - 1 



dk- 



(l + fc2)2 

e 



2Qk^ 



(6.10) 



Both the Tamm-Dancoff approximation (|6.5p and the 
calculations with Eq. (|6.10|) are shown in Fig. [S] Our 
curve is very close to Larsen's one although they lead to 
different values of the edge point momentum {Qc""^ ~ 
1.25). On the other hand, Larsen's curve tends to a 
higher (and wrong) value of the polaron energy at the 
edge point and this is the main source of the discrepancy 
in calculating Qc- However, if we recall that the correct 
edge point value 1 + E{a, 0) is the rigorous upper bound 
for the dispersion (see Eq. (jl.ip ). we have to find the 
momentum at which Larsen's curve reaches it. This way 
we found practically the same Qc ~ 1.20 as given by our 
fitting formula. Thus, we arrive at the same results for 
Qc although Larsen's Eq. (|6.10p and our Eq. (|6.6p do not 
look similar and were obtained in different ways. This 
gives hope that the value of Qc — 1.20 foimd for a = 1 
fits the exact one quite well. 



o 



- Upper bound z(a,Q}, 1 phonon 

- Upper bound z(a,Q}, 2 phonons 
Tamm-Dancoff f^Ja,Q} 

- Larsen fJa,Q) 

- Our fitting curve f(a,Q) 
l\1onte-Carlo 
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FIG. 8: Dispersion of 3D polaron for a = 1 



VII. CONCLUSIONS 

We have presented here the upper and lower bounds 
for the polaron dispersion in D = 1, 2, 3. At small values 
of the electron-phonon coupling constant a the gap be- 
tween our estimates is so narrow that we report in reality 
numerically exact results. 

For intermediate values of a w 1 we proposed the fit- 
ting formula which is simple to use in numerical calcu- 
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lations. Its validity is proved by comparing with the di- 
rect three-phonon calculations for the IZJ-polaron. For 
D = 2, 3 the comparison with the previously obtained re- 
sults is also made, which allows us to conclude that the 
proposed formula fits the exact dispersion quite well in 
all space dimensions. 
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